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We calculate the decay rate of the hghtest heavy Majorana neutrino in a thermal bath 
^ using finite temperature cutting rules and effective Green's functions according to the 

hard thermal loop resummation technique. Compared to the usuai approach where ther- 
mal masses are inserted into the kinematics of final states, we find that deviations arise 
through two difTerent leptonic dispersion relations. The decay rate difi'ers from the usuai 
approach by more than one order of magnitude in the temperature range which is in- 
, teresting for the weak washout regime. This work summarizes the results of Ref. [T] to 



which we refer the interested reader 
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1. Introduction 



' Leptogenesi^ is an extremely successful theory in explaining the baryon asymmetry 

of the universe by adding three heavy right-handed neutrinos Ni to the standard 
model, 

5C = md^-f'^N, - A,,^aAf,0^£a - ^M.iV.iV^ + h.c. , (1) 

with masses Mi at the scale of grand unified theories (GUTs) and Yukawa couplings 
Xv^ia similar to the other fermions. This also solves the problem of the light neutrino 
masses via the see-saw mechanism without fine-tuning.^ 

The heavy neutrinos decay into lepton and Higgs boson after inflation, the decay 
is out of equilibrium since there are no gauge couplings to the standard model. If 
the CP asymmetry in the Yukawa couplings is large enough, a lepton asymmetry is 
created by the decays which is then partially converted into a baryon asymmetry 
by sphaleron processes. As temperatures are high, interaction rates and the CP 
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asymmetry need to be calculated using thermal field theorjBl rather than vacuum 
quantum field theory. 



2. Hard Thermal Loops and Thermal Masses 

When using bare thermal propagators in TFTJS one can encounter IR singularities 
and gauge dependent results. In order to cure this problem, the hard thermal loop 
(HTL) resummation technique has been inventedP^ If g is the coupling to the 
thermal bath, then for soft momenta K < gT, resummed propagators are used. For 
a scalar field e.g. this reads 

.A* = .A + .A(-.n).A + . . . ^ ^ Q^-Ll-U 

The self energy II gT then acts as a thermal mass m^j^ = Il and gives a correction 



3. Decay and Inverse Decay Rate 

Since we are interested in regimes where both the Higgs boson and the lepton mo- 
mentum can be soft, we resum both propagators. The HTL resummation technique 
has been considcred in Ref. |8]for the case of a Dirac fermion with Yukawa coupling. 

In order to calculate the interaction rate T oi N ^ £<j), we cut the N self energy 
and use the HTL resummation for the lepton and Higgs boson propagators. Since 
Ai/,iQ ^ 1, it is justified to neglect the coupling of the neutrino to the thermal bath. 
According to finite-temperature cutting rulesp^ the interaction rate reads 

r(P) = --^tr[(/' + M)ImS(P)]. (3) 

At finite temperature, the self-energy reads 

, r a3i. 

nP)^-g^T 7^PLS*{K)PnD*{Q), (4) 

feo = i(2n+l)7rT ^ ' 

where Pl and Pr are the projection operators on left- and right-handed states, 
Q = P — K and we have summed over neutrino and lepton spins. We also sum over 
the two components of the doublets, particles and antiparticles and the three lepton 
fiavors, such that g^ ~ 4(Aj,Aiy)ii. 

The HTL-resummed Higgs boson propagator is D*{Q) = 1/{Q^ — ni^), where 
m^/T^ = (3/16 5I + 1/16 gì- + 1/4 + 1/2 A) is the thermal mass of the Higgs 
boson. The couplings denote the SU(2) coupling 32, the U(l) couphng gy, the top 
Yukawa coupling yt and the Higgs boson self coupling A, where we assume a Higgs 
boson mass of 115 GeV. The other Yukawa couplings can be neglected since they 
are much smaller than unity and the remaining couplings are renormalized at the 
first Matsubara mode 2ttT as explained in Ref. IH 
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The effective lepton propagator in the hehcity-eigenstate representation is given 



where 



S*iK) = -A+(X)(7o - k • 7) + -A_(X)(7o + k • 7), 



-ko ± k - 



~k 



±1 



2fc 



k ko 
— In ■ 



ko — k 



(5) 



(6) 



and mj/T^ = (3/32 g| + 1/32 g^,). 
The trace can be evaluated as 

tT[{fi + M)PLS*{K)PH]= A+{po^PV) + A-{po+pv), (7) 

where rj — p ■ h/pk is the angle between neutrino and lepton. We evaluate the sum 
over Matsubara frequencies by using the Saclay methodf^ For the Higgs boson 
propagator, the Saclay representation reads 



D*{Q) 



1 



(8) 



/o ^"^g 

where fì~\/T, 713(0;^) = l/(e'^''^ — 1) is the Bose- Einstein distribution and 
(jj"^ — + m?^. For the lepton propagator it is convenient to use the spcctral repre- 
sentatiorP^l 



dr'e 



'„koT' 



duj p±{uJ, k)[l — np{Lo)\e 



(9) 



where ni?(a;) — l/(e"" + l) is the Fermi-Dirac distribution and p± the spectral 
density!^ 

The lepton propagator in Eq. ([5]) has two difFerent poles for 1/A± ~ 0, which 
correspond to two leptonic quasiparticles with a positive (A+) or negative (A_) 
ratio of helicity over chirality.^*"^® The spectral density p± has a contribution from 
the poles and a discontinuous part. We are interested in the pole contribution 

pt^'i^^ k) = ^^^^^(^('^ - '^i) + '^(^ + ^t)), (10) 

where uj± are the dispersion relations for the two quasiparticles, i.e. the solutions for 
1/A ± (w±, k) = 0, shown in Fig. [T] (a). An analytical solution for Ld± can be found 
in the appcndix of Ref. [H One can assign a momentum-dependent thermal mass 
'm,±{kY = i^±(fc)^ — fc^ to the two modes as shown in Fig. [T] (b) and for very large 
momenta the heavy mode m+ approaches Wf, while the light mode becomes 
massless. 

After evaluating the sum over /co, carrying out the integrations over r and r' 
and integrating over the pole part of p± in Eq. (fTO| . we get 



T^D*Ai 

ho 
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2rai 



l + riB — np 
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(11) 
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where ub = ns^i-^q) and np = np{uj±) or ni?(cj^), respectively. 

The four terms in Eq. pT|) correspond to the processes with the energy relations 
indicated in the denominator, i.e. the decay N — > </)£, the production N<j) — ^ £, the 
production N£ — ?> <j> and the production of Nicf) from the vacuum, as well as the 
four inverse reactionsP We are only interested in the process N o (j)£, where the 
decay and inverse decay are illustrated by the statistica! factors 1 + ub — np = 
(1 + rì_B)(l — np) + nsnp, given by the first terra of Eq. (jlip . 

For carrying out the integration over the angle 77, we use 

1 



Im- 



Po - 

After integrating over 77 we get 
1 



•q Kp 



(12) 



T{P) 



167rpoP 



E 



dfc— \M±{P,K)Y[l + nB{ujq±)-np{Lo±)] 

-1<?7±<1 ^± 



1 



(13) 



dfc Aq {2nfS%P -K^Q) \M±{P, K)Y[l + ub - np], 



where ujq± = po — w±, we only integrate over regions with — 1 < ?7 < 1, dfc 
d''fc/((27r)^2 fco) and dq analogously and the matrix elements are 

.,2 fc2 



\M±{P,K)\' =g' 



2mi 



-uj± {pqTPV±) 



(14) 



In order to compare our result to the conventional approximation,^ we do the 
same calculation for an approximated lepton propagator S*ppj.^^{K) = l/(-^ — nii). 



//-2 



This amounts to setting uj"^ — kF' + mj, u}q = pa — uj and we get = ^{M 

nij — as matrix element. 

This result resembles the zero temperature result with zero temperature masses 
m£, m0. The missing factor 1 + nB — np = (1 + nB)(l — np) + riBnp accounts for 
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(b) 

Fig. 1. (a) The two leptonic dispersion relations compared with the standard dispersion relation, 
(b) The momentum-dependent quasiparticle masses rrii^ = upj^ — k^. 
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the statistical distribution of the initial or final particlcs. As pointed out in more 
detail in Ref. ITj we have shown that the approach to treat thermal niasses hke zero 
temperature masses in the final state'* is justified since it equals the HTL treatment 
with an approximate lepton propagator. However this approach does net equal the 
full HTL result. 

4. Decay Density 

The quantity which enters the Boltzmann equations is the decay density integrated 
over ali neutrino momenta. In equilibrium it reads 

- J (S^^-(^)^- - i /Jdi.i.,/-r,, (15) 

where E — pq, f'^{E) = [exp{E/3) — 1]^^ is the equilibrium distribution of the 
neutrinos and Ti, = [1 - f^^iE)] T. 

In Fig.[5J we compare our result to the conventional approximationP We evalu- 
ate the decay rate for Mi — 10^° GeV and mi — (A^Aj^)ii v^/Mi = 0.06 eV, where 
V — 174 GeV is the vacuum expectation value of the Higgs field. 




z = T/M 

Fig. 2. The neutrino decay density with the one lepton mode approach 70 and the two-mode 
treatment 7± for Mi = 10^" GeV and mi = 0.06 eV. The thresholds for the two modes (+), (-) 
and one mode (0) are indicated. 

In the one-mode approach, the decay is forbidden when M < mg + m^. Con- 
sidering two modes, the phase space is reduced for the positive mode due to the 
larger quasi- mass and at M = to+(c!o) -I-to^, the decay is only possible into leptons 
with small momenta, thus the rate drops dramatically. The decay into the nega- 
tive, quasi-massless mode is suppressed due to its much smaller residue. However, 
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the decay is possible up to M = m^. These rates difFer from the one mode ap- 
proach by more than one order of magnitude in the interesting temperature regime 
of z = T/M > 1. 



5. Conclusions 

As discussed in detail in Ref. [ITj we have, by employing HTL resummation and 
finite temperature cutting rules, confirmed that treating thermal masses as kine- 
matic masses as in Ref. |4] is a reasonable approximation. We have calculated the 
decay density of the hghtest heavy Majorana neutrino and its behavior can be ex- 
plained by considering the dispersion relations uj± of the Icpton modes and assigning 
momentum-dependent quasi-masses to them. The thresholds for neutrino decay re- 
ported in Ref. [Jj are shifted and the decay density shows deviations of more than 
an order of magnitude in the interesting temperature regime T/M ^ 1. In order to 
arrive at a minimal consistent treatment, also the decay (j) — > N£ at high temper- 
atures needs to be included as weh as Higgs boson and neutrino CP asymmetries 
which are corrected for lepton modes. 

This contribution summarizes the results of an earher worl^il and we refer the 
interested reader to the more elaborate treatment there. 
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